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1 Introduction 



The high precision of the experimental measurements at the LHC calls for an evaluation 
of the jet, heavy-quark, vector-boson and Higgs-boson production rates at hadron colliders 
which is at least as precise. In recent years, we have witnessed fast progress in the evaluation 
of the production rates mentioned above, in association with many jets, at next-to-leading 
order (NLO) accuracy in the strong coupling constant Og. That rested upon efficient 
methods to compute one-loop amplitudes with many legs, and upon subtraction algorithms 
to evaluate QCD cross sections at NLO at the partonic level [1-5], or at the hadronic 
level, through interfaces, like MC@NLO [6] or POWHEG [7-9] with parton-shower event 
generators. Such algorithms are based on the universality — i.e., on the independence from 
a specific scattering process — of the infrared emissions. 

Lately, the evaluation of production rates at next-to-next-to-leading order (NNLO) 
accuracy has received a lot of attention. Fully differential cross sections for vector-boson [10, 
11], Higgs-boson [12, 13], diphoton [14] and Higgs-vector-boson [15] production have been 
evaluated, and the computation of the NNLO corrections to top-pair production is well 
under way [16-18]. Much work has gone into the lay-out of a general subtraction algorithm 
to compute cross sections at NNLO accuracy [5, 19-37]. The antenna scheme [21] has 
yielded the evaluation of total rates [38-41] and event shapes [42-46] in electron-positron 
annihilation, leading to precise determinations of the strong coupling [47-50]. 

No subtraction algorithm to compute cross sections at hadron colliders at NNLO ac- 
curacy, though, has been devised yet. In order to do that, one must define subtraction 
terms that properly regularise the real-emission phase-space integrals and then one must 
combine the integrated form of those counterterms with the virtual contributions, so as to 
cancel the infrared divergences of the loop amplitudes, in such a way that the cancellation 
of both the kinematic singularities in the real-emission pieces and the explicit e-poles in 
the virtual pieces be local. This implies that the subtraction terms and the real-emission 
contributions must tend to the same value in all kinematic limits where the latter diverge, 
and that the cancellation of explicit e-poles between the integrated subtraction terms and 
the virtual contributions must take place point-wise in phase space. In particular, that im- 
plies that it is possible to write the integrated counterterms in such a way that they can be 
explicitly combined with virtual contributions, before phase-space integration. Practically, 
the locality of the subtraction scheme is also important to ensure good numerical efficiency 
of the algorithm. In broad outline, we remind how this occurs at NLO. After fixing the 
leading-order m-jet cross section as the integral of the fully differential Born cross section 
do"^ of m final-state patrons over the available m-parton phase space defined by the jet 
function J^, 





the NLO correction. 
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can be rewritten as a sum of finite integrals, 



a 



NLO 



m+1 



^NLO 

'm+1 



NLO 
'm ) 



where each of the two terms on the right-hand side 

da^NLO 
do", 



'm+1 
NLO 



^^m+l^m+l — do'j^^lJm 



'm 



R,Ai 

m+1 I '^m 



e=0 



e=0 



(1.3) 

(1.4) 
(1.5) 



is integrable in four dimensions by construction [1-5]. 

Likewise, we can write the NNLO correction to the cross section as a sum of three 
contributions, the doubly unresolved, the one-loop singly unresolved, and the two-loop 
doubly virtual terms. 



NNLO 



m+2 



'^'^ m+2'' m+2 



+ [ da^+^Jrn+i+ [ dal'^Jm, (1.6) 
J m+1 Jm 



and rearrange it as follows. 
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where each of the three terms on the right-hand side. 
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(1.7) 

(1.8) 
, (1.9) 

}Jm 1 

(1.10) 



is integrable in four dimensions by construction [20, 22, 23]. 

The counterterms which contribute to dcj^^2*^ and to dcr^^^j'^ were introduced in 
Refs. [22] and [23], respectively. The integral of the real- virtual counterterms (the last 
two terms of Eq. (1.10)) was performed in Refs. [25-27]. The integral of the iterated 
singly unresolved counterterm (the third term of Eq. (1.10)) was performed in Ref. [31]. 
In this paper, we compute the integral of the collinear-type contributions to the doubly 
unresolved counterterm (the second term of Eq. (1.10)). The soft-type contributions to the 
same integral are presented in a companion paper [51]. 



2 Notation 



In this paper, we use the notation introduced in Ref. [31], which we recall in this Section. 
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2.1 Matrix elements 

We consider processes with coloured particles (partons) in the final state, while the initial- 
state particles are colourless (typically electron-positron annihilation into hadrons). Any 
number of additional non-coloured final-state particles are allowed, too, but they will be 
suppressed in the notation. Resolved partons in the final state are labeled with letters 
chosen from the middle of the alphabet, i, j, k, I, . . . , while letters r, s denote unresolved 
final-state partons. 

We adopt the colour- and spin-state notation of Ref. [2]. In that notation the ampli- 
tude \M) for a scattering process is an abstract vector in colour and spin space, and its 
normalisation is fixed such that the squared amplitude summed over colours and spins is 

IMl"^ = {M\\M) . (2.1) 

In this paper we only need \M) in the Born approximation, denoted by 

Using the colour-state notation, we define the two-parton colour-correlated squared 
tree amplitudes as 

= = (>i(°^|r,rr,,|>[W) (2.2) 

and similarly the four-parton colour-correlated squared tree amplitudes, 

\^tkw/ = ® {T^■Tk,TJ.Tl} ^ (XW|{T,.Tfc,T,.Ta|>lW) , (2.3) 

where we also introduced the Cg) product notation to indicate the insertion of colour 
charge operators between (A^^'^^j and The colour-charge algebra for the product 

ZniTinTkr = T,-Tk is: 

T,-Tk = Tk-Ti if i^k; Tj = Cf^. (2.4) 

Here C/. is the quadratic Casimir operator in the representation of particle i and we have 
C7g = Cf = TR(iV2 - l)/iVc = (A^-^ - l)/(2iVc) in the fundamental and Cg = CA = 2 TrN^ = 
Nc in the adjoint representation, with Tr = 1/2. 

We also use squared colour charges with multiple indices, such as T?, = C/.^ and 
Tf^g = Cf.^^. In such cases the multiple index denotes a single parton with flavour obtained 
using the flavour summation rules: odd/even number of quarks plus any number of gluons 
gives a quark/gluon, or explicitly for the relevant cases at NNLO, 

q + 9 = q, q + q = 9, 9 + 9 = 9, 

q + 9 + 9 = q, q + q + Q = q, 9 + q + Q = 9, 9 + 9 + 9 = 9- (2-5) 

2.2 Cross sections 

In this paper we shall need to use only tree-level n-parton production cross sections, with 
n = m, the Born cross section, and n = m + 2, the so-called doubly real correction. We 
have 
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where N includes all QCD-independent factors and d(j)n{{p}) is the d-dimensional phase 
space for n outgoing particles with momenta {p} = {pi, . . . and total incoming mo- 
mentum Q, 

dMPu ...,Pn;Q) = ll S+ip^) {27,^5^'^ iQ-Y^pA. (2.7) 

i=i ^ ^ t=i ^ 

The symbol Yli{n} denotes summation over different subprocesses and S'{„} is the Bose 
symmetry factor for identical particles in the final state. Then the Born cross section and 
the doubly real correction are simply 

Aaliip}) = dat\{p}) and da^^^lW) = '^^tUiiP}) ■ (2-8) 

The final result will also contain the phase-space factor due to the integral over the 
{d — 3)-dimensional solid angle, which is included in the definition of the running coupling 
in the MS renormalisation scheme^, 

^-y(2^-f(r^- ^'-'^ 

2.3 Momentum mappings and phase-space factorisation 

The subtraction terms are written in terms of the momenta obtained by the mappings of 
Ref. [22], 

{p}m+2 ^ WL^^^ , (2.10) 

where — > may label a triple collinear, double collinear or soft collinear mapping^. In 
particular, the soft collinear mapping is obtained by the iterated application of a basic 
collinear mapping, followed by a basic soft mapping, 

{p}m+2 ^ MSi ^ {P}^^''^ , (2.11) 

while the former two cases do not have such an iterated form. As the above notation sug- 
gests, the final set of m momenta are denoted by tildes, while hats indicate the intermediate 
set of m -|- 1 momenta. In kinematic expressions where only the label of a momentum is 
displayed (we shall discuss several examples below), the tilde is inherited by the label, and 
we write for instance i , ir and irs , where the last two label single momenta. However, 
since these mappings affect only the momenta, but not the colour and flavour (apart from 
the flavour summation rules of Eq. (2.5)), we shall omit the tilde from flavour and colour 
indices. 



""^In the MS renormalisation scheme as often employed in the literature, the definition of the running 
coupling includes the factor Se = {4nye~'^''^ . In a computation at NLO accuracy, the two definitions lead 
to the same expressions. At NNLO they lead to slightly different bookkeeping of the IR and UV poles at 
intermediate steps of the computation, our definition leading to somewhat simpler expressions. Of course 
the physical cross section does not depend on these details. 

^Integration of the subtraction terms obtained with the double soft mapping are presented in the com- 
panion paper [51]. 
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We also use labels such as (ir) to denote a single momentum that is simply the sum 
of two momenta, p^.^^ = pf + Pr- 

The momentum mappings are chosen as to lead to an exact factorisation of the phase 
space. For the triple collinear phase-space mapping, we have, 

d<Pm+2{{p}m+2;Q) = d0™ , Q) [dj^;^ (pr , Ps , P^rs^Q)] , (2-12) 

for the double collinear mapping, 

d</.^4-2(M^+2; Q) = d^„rmt''\Q) [dp^S'^^iPr^Ps, P„. P js', Q)] , (2-13) 

while for the iterated collinear and soft mapping appearing in Eq. (2.11), 

d<A^+2(W^+2; Q) = dcPrnmi^''^^) [dpg+, (p, , p,, ; Q)] [dp^liPs; Q)] ■ (2.14) 
2.4 Kinematic variables 

The following types of kinematic variables are used to write the doubly unresolved sub- 
traction terms. 

• Two-particle invariants, such as 

Sir = 2pi-pr or SiQ = 2pi-Q. (2.15) 
Two-particle invariants scaled with are denoted by yij = Sij/Q'^. 

• Momentum fractions Zi^r for the splittings pir Pi + Pr and zi^j-s for the splittings 

Pirs Pi +Pr+Ps, 

Zir = and Zira = , (2-16) 

ViQ + VrQ ' ViQ + VtQ + VsQ 

with Zi^r + Zr^i = 1 and Zi^rs + Zr^si + Zg^ir = 1 {zr,si and Zs^ir are computed by obvious 
permutations). 

• We also use the eikonal factors, 

S,u{s) = ^. (2.17) 

As mentioned above, the sum of two momenta is often abbreviated with the two indices 
in parentheses e.g., +Pr = P^i^)-: which is also used in other occurrences, such as 

S(ir)k = Sik + Srk , <S(^ir)k{s) = J — 7 • (2-18) 
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Finally, we express the integrated subtraction terms as functions of the following (com- 
binations of) invariants: 



x~ = 



and 



Y~-r 



yiqy jc 



(2.19) 



In the rest frame of Q^, y^g = 2Ei/Q is simply the scaled energy of parton i, while g 
is related to Xij-, tlie angle between momenta and p^, by cosxij = 1 — q- 



3 Integrating the doubly unresolved approximate cross section 



3.1 The integrated approximate cross section and insertion operator 

The doubly unresolved approximate cross section times the jet function is defined as, 

da^^^ QJ^=M ^^,^^,i^{p})-l—A^^)\Mt\,{p)\' j„(p) , (3.1) 



{m+2} 



'{m+2} 



where the notation reminds us that the set of momenta entering Jm is different from 
term to term, 

I-^212(P)I'0^™(P) = 
1 

6 



EE E 



(s,jr) ~\ 



(3.2) 



All terms on the right-hand side of Eq. (3.2) are defined in Ref. [22]. 

Using that all momentum mappings lead to the exact factorisation of phase space, 
d^m+2({p}; <3) = d(^m({p}; (5)[dp2;m]; and that the jet function does not depend on 
the variables of the factorised two-parton measure, [dp2;m] > we can compute the integral of 
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Eq. (3.2) over the two-parton factorised phase space independently of J, 




\ j^i,r,s k^j 





i,k \ j,l 



(3.3) 



In order to make the integrated subtraction terms dimensionless in color space, we have 
factored out quadratic Casimir operators in Eq. (3.3). 

Eq. (3.3) is not yet in the form of an m-parton contribution times a factor. In order to 
obtain such a form, we must still perform summation over unresolved flavours, rewriting the 
symmetry factor of an m + 2-parton configuration to the symmetry factor of an m-parton 
configuration. The counting is performed as in the Appendix of Ref. [31]. As a result of 
the summation over the unresolved flavours, we obtain functions — the flavour-summed 
integrated counterterms — denoted by (X^'^^)^"'''^'", which are specific sums of the non 

flavour-summed integrated subtraction terms, [X^^^]^ji'''^"' . Schematically, we may write, 



The flavour indices on the left- and right-hand sides of this equation need not match: 
the non flavour-summed functions on the right-hand side carry dependence on unresolved 
flavours, while the flavour-summed functions on the left do not, by definition. 

After summing over unobserved flavours, the integrated doubly unresolved approxi- 
mate cross section can be written as. 



where the insertion operator (in colour space) has three contributions according to the 




(3.4) 




(3.5) 



8 



possible colour structures, 



27r 



ii. 



2 \ e 



E 



+ EE 



j(o),(i,0 



(o),(i,/) ^2 



T,Ti (3.6) 



i,j,k,l 



{TiTk,TjTi}} , 



with = C/^ {Cg = Cf, Cg = Ca) as in Eq. (2.4). 

In terms of flavour-summed integrated counterterms discussed above, the kinematic 
functions of the insertion operator can be written as, 



p(0) 



c 



(0) 
2,ij 



cs 



(o),(i,0 

2,1 



(0) 



C, 



CS 



(0) ■ 

ir;js 

(0) 

ir;s 



Ji 



S(o) 
S(o) 



i 



p q(o) I -I- ( f rs s^*') 

^irs^rs I r \ irs^-^ir;s^rs 
' Ji 



+ (C- . 



fi 



Uf3 



(3.7) 



On the right-hand side of Eq. (3.7), the flavour-summed functions depend on the kinematics 
through variables of the type xi and Yij^Q- The latter dependence stems from integrating 
an eikonal factor which is always multiplied by a colour-connected squared matrix element. 
In order to make the results more transparent, we hid the arguments of the functions, but 
kept the relation to the colour-connected matrix elements, shown as upper indices. 

3.2 Flavour-summed integrated counterterms 

On the right-hand side of Eq. (3.7), the flavour-summed functions are the following combi- 
nations of the integrated subtraction terms. 



1. Triple collinear: 



C 



(0) 



(0) 



a 



(3.8) 



irsiaqq 



The factor (nf — 1), which appears after summing over the unresolved flavours, can 
be traded for nf by introducing the integrated subtraction term [C^yM'^J, defined by 



n 



(0)1 



irsiqqq 



(0)1 (id) 



irsiqqq 



(3.9) 
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This definition matches the decomposition of the same- flavour triple splitting function 
computed in Ref. [52]. Thus, we obtain 

(Cl°l)^ = ^[C!2],,, + nf [C\%,,,, + l[Ctl]'^;i^ , q'^q. (3.10) 

2. Triple collinear - soft collinear: 

f io)\ 1 (0) (0) (^-^^^ 

3. Double collinear: 

(C^O) \ _l[n(0) 1 

/p{0) ^ _ lrp{0) n !^rp(0) 1 

\yir;jsj^^~^v^ir;js\gg;Qg^ 2 ^ ir-jskmg ^ V-J-J-^J 
\^ir;jsj^^ ~ gl^ir;j.s\gg;gg ^ 2 ^ ir;js\qq;q q 

4. Double collinear - soft collinear: 

(o)\ 1 (0) (0) ^^--^^^ 

y^ir;jf^ir\s ) gj- ~ 2 f^*'"'^^^*''''^-^^^ i^ir;jl~^ir;s\qq ' 

i.e., it is independent of the flavour /. 

5. Soft collinear: 

fcs(o)y^'')^rcs(o)i(i,0 

\^^ir;sl l^^ir;s\qg ' 

r-n 1 (3-14) 

/pq{0) V^''-" _ ^ rpq;(0) nOV) rpc;{0) 
[^ir;s)^ - 2 ^^ir;Jgg + i^ir;Jgg • 

6. Triple collinear - double soft: 

^Cjj,^S^^^^ ^ = 2\-^irs^r^s]fgg + '^f [CirsS^s^J/gg . (3.15) 

7. Triple collinear - soft collinear - double soft: 

(c,„cs,„s<»))^ = [c.„cs„,sS)i. (3.e) 

i.e., it is independent of the flavour /. 
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8. Double collinear - double soft: 

\^ir;js^rs) . . 77 [^jry's^^^''] , (3.17) 
^ ^ /1/2 ^ 

i.e., it is independent of the flavours fi and f2- 

9. Double soft: 

kj^a I r, I'-'rs Joo ^ t L'-'rs loo 



(3.18) 



10. Soft collinear - double soft: 

(3.19) 

i.e., it is independent of the flavour /. 

In this paper, we compute the collinear-type counterterms, which do not involve S^J , 
i.e., those in Eqs. (3.10)-(3.14). The integration of the subtraction terms which are obtained 
by the double soft mapping is presented in a companion paper [51]. 

4 Integrated counterterms 

In this section we define the integrated counterterms and compute them in terms of master 
integrals. 

4.1 Integrated triple collinear counterterm 

The integral of the triple collinear counterterm is 

V / J2 Sjrs [J-irs) (4.1) 

^ '^/i/r/s ^^,isi ^s,iri Sir, Sis, Srs, c) f ip^Q, Oiirs, dijTL, 6)) . 

Here Pf^j^f^ {zi,rs, Zr,is, Zs,ir, Sir, Sis, Srs] e) are the complete^ spin-averaged three-parton split- 
ting functions, recalled in Appendix B, where we explain why some of these are different 
from those presented in Ref. [52]. For consistency of the complete subtraction scheme one 
must use the forms given in Appendix B. The function f(a(),airs,d{m,e)) is defined and 
its role is explained in Appendix A. The other subtraction terms discussed in this paper 
will also contain such harmless modifications as compared to the original ones in Ref. [22]. 

Following the decomposition of the triple collinear splitting functions into abelian and 
non-abelian pieces [52], we also decompose the integrated functions [C-|?^]qgg and [C^|?]]ggg 
likewise, 

rp(0)l _ rp(0)n(ab) , rp(0)n(nab) , rr.(0)n _ _ rp(0)n (ab) rp(0)n (nab) . 

V^irs\Q99 ~ l^irslqgg ' V^irsiqgg ^liu [^irs\gm ~ l^irslgqg ' l^irslgqg ■ l^-^J 



As opposed to abelian and non-abelian parts, see Appendix B.2. 
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Table 1. Non-zero coefficients of the Xj^'^j'''™' (x, e; ao, do) integrals. 
They are computed in Appendix C. The result can be written as 

[c!2]/.^^ = «^^/. E E -i°}3?:r:^a'''"^fe,^;«o,do), (4.3) 

n=l j,k,l,m 

where the constants CLf-f^j^ denote colour-factor ratios as follows: 

_ Tr (id) 
dqg'g' — ! CLqqq 

(ab) C'fTr (nab) 
n - = n - 

The non-zero coefficients are listed in Tables 1-5, while the integrals X2^''^'''"''(x, e; (^o) 
(n = 1, . . . , 5) are defined and computed in Appendix C. 
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-12 

6 



Table 2. Non-zero coefficients of tlie 22c'2^'"^\^t ^'-i c^Oi do) integrals. 



4.2 Integrated double collinear counterterm 

The integral of the double collinear counterterm is 



1 1 



SirSjs T'^^T'jg 



X Pfj^izr,i; e)Pyjj^{zs,j\e) f{ao, a^. + ajs, d{m, e)) . 



(4.5) 



In Eq. (4.5), p!^^j^{zr^i; e) and PjPj^{zr^i; e) are the spin-averaged Altarelli-Parisi splitting 
functions in d dimensions, as recalled in Appendix B. 

The integrated counterterm is computed in Appendix D. The result is 

2 

k,l=-l 



rp(0) 1 -n n \^ (0),fc (0),/^(fc,0. ^ ^ 



(4.6) 



where the constants a/^/^ denote the colour factor ratios 



''qq 



Ik 

Ca 



0.qn — 0,„„ — 1 , 



, U.qg — Ugg 



(4.7) 



and the coefficients Cf^f^ are those already defined in Ref. [27] stripped of their colour 
factors, as presented in Table 6. 
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j 


k 


/ 


m 




^(0),(id) 
qqq;3 


(0),(ab) 


JO), (nab) 


g(0)^(ab) 


JO), (nab) 
911 t-J 







-1 





-1 


2 








-8 





2 


-24 





-1 


1 


-1 


-2 


-4 + 4e 





2 - 2e 





-2 


6-66 





-1 


2 


-1 


2 








-2 + 2e 





8 

1 - e 


-24 


-1 








-1 





2-2e 





-4 





1 


-12 


1 

— i 


u 


1 
i 


1 

— 1 


n 
u 


■4:6 


n 
u 


O — 6 


u 


1 - e 


1 


1 

— - i 


u 


9 


1 

— 1 


n 
u 


Z — Zfc 


n 
u 




u 


4 

1 - e 


1 9 

— IZ 


1 

— i 


u 


Q 
O 


1 

— 1 


n 
u 


n 
u 


n 
u 


n 
u 


n 
u 


2 

1 - e 


D 


u 


9 


1 
i 


1 

— 1 


A 


n 
u 


n 
u 




n 
u 




1 9 _l_ 1 9,; 
— IZ -p IZfc 


1 


-2 


1 


-1 


-4 








4 - 4e 





-4 


12 - 12e 


-1 


-1 





-1 





-4 





2 











_x 


_x 


I 


_1 





—2 — 2f 





_2 







— 12 


-1 


-1 


2 


-1 





-2 + 2e 





1 - e 





2 

1 - e 


6 


-1 


-1 


3 


-1 

















2 

1 - e 


-6 





-2 


2 


-2 


-2 








2 - 2e 





-2 


6-6e 



Table 3. Non-zero coefficients of tlic l2Q'^'"^\x,t;aQ,dQ) integrals. 



3 


k 


I 


m 


(0) 


^(0),(id) 
^qqq-A 


(0),(ab) 

qgg-A 


(0),(nab) 

^qgg-A 


(0),(ab) 
'-qqq-A 


(0),(nab) 
'-gqq;A 


^ggg-A 


-1 





-1 


-1 








4 











12 





-1 


-1 


-1 








4 


-2 








6 


-1 


-1 


-1 


-1 








4 


-2 








6 



Table 4. Non-zero coefficients of the Ijc 4' ™^ (2^' o^o, c^o) integrals. 



4.3 Integrated soft coUinear counterterms 

There are three types of soft collinear counterterms. 
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j k I m 


(0) 


„(0),(id) 

qqq;5 


(0),(ab) 

qgg;5 


(0),(nab) 

qgg;5 


(0),(ab) 

gqq;5 


(0),(nab) 

gqq;5 


ggg;5 


-1 -1 -1 -1 





2 














3 


Table 5. Non-zero coefficients of the I2q'^'"^\x 


, e; do) integrals. 







„(o),fc 
'^qg 




^gg 


-1 


2 





4 





-2 


1 


-4 


1 


1 - e 


2 

1 - e 


2 


2 





2 

1 - e 


-2 



Table 6. Values of the cf^j'^ coefficients that appear in Eqs. (4.6), (4.9), (4.f0), (4.12) and (4.14). 



1. Soft collinear: 



[CS, 



(0) 



(47r^2 X 2 



-Q 



2e 



[(^Plyl+l(Pr,Pir; Q)] [dp\".l^{ps] Q)] 



X ^■Sjkis) —^Pfj^{zr,i;e) f{ao,air,d{m,e)) f(,yo,ysQ,d'{m,e)) 

^ Sir -I 



(4.8) 



The factor Sjk{s) is the eikonal factor, given by Eq. (2.17) if j,k / (ir), and by 
Eq. (2.18) ife.g., j = (ir). 

The integrated counterterm is computed in Appendix E. The result is 



l=-l 

(4.9) 

if j, k / (ir), and 



(4.10) 



if e.g., j = (ir), where af.f^ and the coefficients are given in Eq. (4.7) and Table 6 
respectively. 
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2. Triple collinear - soft collinear: 



PfjMr,u e) /(ao, air, d{m, e)) f{yo,ysQ,d'{m, e)) 

(4.11) 



2 1 - Zs^ir _l p(0) 



which is computed in Appendix E. The result is 

2 



l=-l 



3. Double collinear - soft collinear: 

N2 N 2 



ICir^jf^S^slMr = {^Q'') l[dpf!l+liPr,Pir;Q)] [dpfliPs; Q)] 

X — ^ — P^^j^{zr,i;e) f{ao,air,d{m,e)) f{yo,ysQ,d'{m,e)) , 

(4.13) 

which is also computed in Appendix E, as 

2 

[Qr;isCsS;l]/,/r = a/,/, Yl ^2^,5(^7;^' %J,q; "0,^0, do, do) • (4.14) 

Z=-l 

5 Results 

After evaluating the integrals as explained in the Appendices, we obtain the kinematic 
functions defined in Eq. (3.7) as Laurent expansions in e. We compute those expansions 
analytically up to 0(e"^), while the remaining coefficients up to 0(6*^) are calculated nu- 
merically via sector decomposition, see [53] and references therein. In obtaining our results, 
we use do = Dq + die and dQ = D'q + d'^e, see Eq. (A. 9). The parameters Dq, di, D'q and 
d'l, as well as ao and yo are left symbolic throughout the analytic computation and sector 
decomposition. 

To present analytic results for the non-soft counterterms up to 0(e~^), we use the 
following abbreviations, 

1qi'>T'i) = ^, 79(^f) = y - ^'^f and Y.{z,N) = lnz ' (^'^^ 

1=1 

where 'Jqinf) is actually independent of the number of light flavours, but introducing the 
flavour dependence formally makes possible a flavour-independent notation. Furthermore, 
the kinematic functions are dimensionless in colour space, hence these rif-dependent 7j(nf) 
constants are related to the constans 7j often used in the literature [54] by jj — ?• Cj7j(nf). 



1. Triple collinear: 



(5.2) 



2. Triple collinear - soft collinear: 

+ 0(e 



( 



p (~c{0) 



(5.3) 



3. Double collinear: 



C 



(0) 
ir-.js 



( ^ J L 

263 



2(lnxi + Inxj) - 7/,("-f) - 7/j (?^f) 



+ 0(e-2) 
4. Double collinear - soft collinear: 



(5.4) 



/-I /-in{0) 



5. Soft collinear: 



CS 



(0) 



(i,0 

I 

h 



+ 0(6-2) _ (5 5) 

^ + ^ ( In Xi + S(2/o, I^o - 1) ) + ^ ( In Y^i-q - jf^ (nj )) 



+ 0(6-2) 



(5.6) 



for j,l ^ i, and 



(0) 



(2;j,yi/.Q) = - 



1 2 

+ ^( Inxj + T.{yo,D'Q - 1] 



g4 g3 



+ ^(lny./;Q-^7/.(nf))+ 0(6-2) 



(5.7) 



for e.g., j = i. 



The remaining coefficients in the Laurent expansion are computed numerically. By 
way of illustration, we present results for the flavour-summed counterterms in Tables 7-19 
for two kinematic points: one relevant for 2-jet production and the other corresponding to 
the fully symmetric configuration of final state momenta in 3-jet production. In terms of 
the invariants introduced in Eq. (2.19), these two phase-space points correspond to 



2-jet : xi = X2 = 1 



and Y-i 



12,1 



1 



2 3 
3-jet symmetric : xi = X2 = X3 = - and Yi2,q = Yis.Q = >23,Q = ^ 



(5.8) 
(5.9) 
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In the numerical computations, we chose the following values for the phase-space cut pa- 
rameters: ao = I, Do = I, Dq = Dq = 3 and di = d[ = —3. In Tables 7-19, a displayed 
error estimate of ±0.000 implies that the numerical uncertainty for that particular value 
is smaller than 10"'^. 



colour 


X 




1/ei 


finite 


1 


1.000 
0.667 


-1.870 ± 0.003 
2.371 ± 0.004 


-22.187 ± 0.020 
-15.435 ± 0.023 


-76.907 ± 0.155 
-93.639 ± 0.156 


Ca 
Cf 


1.000 
0.667 


2.632 ± 0.002 
6.080 ± 0.002 


4.595 ± 0.013 
21.854 ± 0.014 


13.916 ± 0.095 
84.489 ± 0.100 


IfTR 

Cf 


1.000 
0.667 


-0.611 ± 0.001 
-0.881 ± 0.001 


1.159 ± 0.005 
-1.016 ± 0.006 


16.629 ± 0.037 
7.092 ± 0.046 



Table 7. CocfRcicnts of the Laurent expansion of (^C^,,^j (x) for do = 3 — 3e and ao = 1. 



colour 


X 




l/.i 


finite 


1 


1.000 
0.667 


3.664 ± 0.003 
11.993 ± 0.004 


-8.987 ± 0.029 
19.817 ± 0.032 


-42.269 ± 0.207 
30.542 ± 0.224 


Ca 


1.000 
0.667 


-2.614 ± 0.005 
-3.424 ± 0.006 


-5.230 ± 0.047 
-11.873 ± 0.049 


-3.361 ± 0.345 
-32.821 ± 0.372 


CFrifTa 
Ca 


1.000 
0.667 


-4.798 ± 0.000 
-5.540 ± 0.000 


-17.623 ± 0.003 
-24.021 ± 0.004 


-45.473 ± 0.015 
-74.365 ± 0.020 



Table 8. CoefRcients of the Laurent expansion of ( C^^^ ) (x) for do = 3 — 3e and ao = 1- 



colour 


X 


Y 






finite 


1 


1.000 
0.667 


1.000 
0.750 


-7.311 ± 0.003 
-13.938 ± 0.004 


-3.228 ± 0.009 
-20.393 ± 0.019 


-9.300 ± 0.027 
-27.573 ± 0.066 



Table 9. Coeffieients of the Laurent expansion of |^CS,-,..gj i^^Y) for do ^ (I'q = 3 — 3e and 
ao =ya = 1- 
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rolnnv 


X 


Y 






finitG 


1 


1.000 
0.667 


1.000 
0.750 


-9.205 ± 0.003 
-16.113 ± 0.004 


-8.103 ± 0.009 
-27.124 ± 0.019 


-16.717 ± 0.028 
-40.577 ± 0.067 


Ca 


1.000 
0.667 


1.000 
0.750 


4.454 ± 0.001 
5.017 ± 0.001 


14.204 ± 0.004 
18.479 ± 0.008 


29.039 ± 0.020 
44.490 ± 0.046 



Table 10. Coefficients of the Laurent expansion of (CS^yi* 
"0 = yo = 1- 



{x, Y) for do = c?o = 3 — 3e and 



colour 


X 


Y 


1/.2 


l/ei 


finite 


1 


1.000 
0.667 


1.000 
0.750 


-4.874 ± 0.002 
-10.330 ± 0.003 


2.910 ± 0.007 
-8.326 ± 0.018 


5.559 ± 0.024 
9.292 ± 0.087 



Table 11. Coefficients of the Laurent expansion of 




= 3 — 3e and 



"0 = yo = 1- 



colour 


X 


Y 


1/.2 




finite 


1 


1.000 
0.667 


1.000 
0.750 


-6.126 ± 0.002 
-11.839 ± 0.003 


0.812 ± 0.007 
-12.094 ± 0.018 


7.619 ± 0.025 
6.585 ± 0.089 


Ca 


1.000 
0.667 


1.000 
0.750 


3.005 ± 0.000 
3.517 ± 0.001 


7.201 ± 0.002 
11.052 ± 0.007 


3.080 ± 0.017 
16.471 ± 0.045 



Table 12. Coefficients of the Laurent expansion of ^CS-^l'gj ix,Y) for do ~ ~ 3 ~ 3e and 
"0=2/0 = 1- 



colour 


X 




l/.i 


finite 


1 


1.000 
0.667 


3.437 ± 0.001 
5.886 ± 0.001 


-6.996 ± 0.004 
-6.784 ± 0.005 


-61.629 ± 0.015 
-100.401 ± 0.019 



Table 13. Coefficients of the Laurent expansion of |^Cj^^CS-^.\, j (x) for c?o = c'o = 3 — 3e and 
ao = yo = 1. 
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rolonv 


X 






finite 


1 


1.000 
0.667 


4.048 ± 0.001 
6.633 ± 0.001 


-6.973 ± 0.004 
-5.975 ± 0.005 


-66.612 ± 0.016 
-103.489 ± 0.019 


Ca 


1.000 
0.667 


-1.556 ± 0.000 
-1.826 ± 0.000 


-1.601 ± 0.001 
-3.445 ± 0.001 


8.364 ± 0.011 
2.731 ± 0.011 



Table 14. Coefficients of the Laurent expansion of (^Cj^^CS^^.j, j {x) for do = rfg = 3 — 3e and 
ao = 2/0 = 1- 



colour 


X 


y 




1/^^ 


finite 


1 


1.000 
0.667 


1.000 
0.750 


8.311 ± 0.003 
13.099 ± 0.002 


11.149 ± 0.010 
24.109 ± 0.011 


40.832 ± 0.028 
61.583 ± 0.031 



Table 15. Coefficients of the Laurent expansion of |^Cj^.j,,CS^^^, j {x, Y) for do = d'^ = 3 — 3e and 
"0 = yo = 1- 



colour 


X 


Y 






finite 


1 


1.000 
0.667 


1.000 
0.750 


10.205 ± 0.003 
15.179 ± 0.002 


16.358 ± 0.010 
30.471 ± 0.011 


51.240 ± 0.029 
75.387 ± 0.032 




1.000 
0.667 


1.000 
0.750 


-4.454 ± 0.001 
-4.825 ± 0.001 


-14.871 ± 0.004 
-17.549 ± 0.004 


-35.687 ± 0.020 
-45.158 ± 0.022 



Table 16. Coefficients of the Laurent expansion of C,^. ,CS "4 (x, Y) for do = d'o = 3-3e and 
ao = 1- 



colour 


Xl 


X2 


l/e^ 




finite 


1 


1.000 
0.667 


1.000 
0.667 


0.162 ± 0.000 
4.408 ± 0.001 


-4.108 ± 0.003 
5.732 ± 0.004 


5.189 ± 0.013 
18.494 ± 0.018 



Table 17. Coefficients of the Laurent expansion of ( C,-,,. ^ ) {xi, X2) for do — 3 ~ 3e and = 1- 

V '-^ / qq 



6 Conclusions 

We have computed the integrals over the two-particle factorised phase space of the collinear- 
type contributions to the doubly unresolved counterterm of the NNLO subtraction formal- 
ism, defined in Refs. [22, 23]. We presented those integrals in terms of parametric rep- 
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folmiT" 




To 


l/e^ 




finite 


1 


1.000 
0.667 


1.000 
0.667 


0.859 ± 0.000 
5.333 ± 0.001 


-3.131 ± 0.003 
8.076 ± 0.004 


5.070 ± 0.013 
22.304 ± 0.019 


Ca 


1.000 
0.667 


1.000 
0.667 


-1.727 ± 0.000 
-2.183 ± 0.000 


-3.683 ± 0.002 
-6.871 ± 0.002 


-3.445 ± 0.008 
-14.493 ± 0.011 



Table 18. Coefficients of the Laurent expansion of ( C^^. ) {xi, X2) for rfo = 3 — 3e and ao ~ 1, 

with (cfJjA {X1,X2) = (cfJjA (.T2,Xi). 
\ / an \ •■' / no 



colour 


Xl 


X2 






finite 


1 


1.000 
0.667 


1.000 
0.667 


1.612 ± 0.000 
6.314 ± 0.001 


-1.855 ± 0.003 
10.780 ± 0.004 


5.670 ± 0.013 
27.243 ± 0.019 


Ca 


1.000 
0.667 


1.000 
0.667 


-3.677 ± 0.000 
-4.589 ± 0.000 


-8.780 ± 0.003 
-15.404 ± 0.002 


-11.184 ± 0.011 
-35.158 ± 0.016 




1.000 
0.667 


1.000 
0.667 


0.222 ± 0.000 
0.222 ± 0.000 


1.637 ± 0.000 
1.884 ± 0.000 


5.932 ± 0.003 
8.057 ± 0.004 



Table 19. Coefficients of tlie Laurent expansion of ( ClJj^ j {xi,X2) for do = 3 — 3e and ao = 1- 



resentations that are suitable for evaluation with sector decomposition. After evaluating 
them, we checked the numerical results with the publicly available code SecDec [55], al- 
ways finding agreement within the numerical uncertainty of the integrations. The soft-type 
contributions are presented in a companion paper [51]. 

By these two papers, we complete the integration of the subtraction terms over the 
unresolved phase spaces, and the computation of the finite cross section in Eq. (1.9) becomes 
feasible for electron-positron annihilation into two and three jets. Although the formalism 
is complete, for a higher number of jets some more work is required, because some integrals 
were evaluated specifically for three-jet kinematics. 
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A Modified doubly real subtraction terms 

We outline a simple modification to the NNLO subtraction scheme presented in Refs. [22, 
23]. Parts of these modifications were introduced previously: those relevant to the singly 
unresolved approximate cross section dcr^^'^^ in Eq. (1.8), and to the approximate cross 
sections in Eq. (1.9), were introduced in Ref. [25], while those relevant to the iterated doubly 
unresolved approximate cross section dcr^^'^^^ appearing in Eq. (1.8) were presented in 
Ref. [31]. 

Recall that the doubly unresolved approximate cross section can be written symboli- 
cally as in Eq. (3.1) 

dcrTi}' = d<t>n.[dp2]Af^\M^Z2\' , (A.l) 

where the doubly unresolved approximation is a sum of terms (see Eq. (3.2)). The precise 
definition of these terms involves the momentum mappings discussed in Section 2.3. All 
such mappings lead to an exact factorisation of the m-|-2-particle phase space, symbolically 
written as 

d(^„+2(M;Q) = d(j),ni{p}m;Q)[dp2:m] ■ (A.2) 

The only feature of the factorised phase spaces [dp2;m] that is relevant presently is that 
they carry a dependence on the number of partons, m, of the form 

a„,,)2(— , (A.3) 
a,r - a,-,)2(— , (A.4) 

yrQ-ysQ + yrs)^"^-'^^'-'^-', (A.5) 

The subtraction terms, as originally defined in Ref. [22] do not depend on the number of 
hard partons, thus the m-dependence of the factorised phase space measures is carried over 
to the integrated counterterms, where this dependence enters in a rather cumbersome way 
(see e.g., Eqs. (A.9) and (A. 10) of Ref. [4]). 

Thus, as in Refs. [25, 31], we reshuffle the m-dependence of the integrated counterterms 
into the subtraction terms themselves, where it appears in a very straightforward and 
harmless way^, through factors of (1 — a) and/or (1 — y) raised to m-dependent powers. 
We gather the results in Table 20, where together with the subtraction terms, we give the 
momentum mappings used to define the term(s) and the expression which multiplies the 
original counterterm to produce the modified one. The function / in Table 20 is defined as 

f{zo,z,p) = @{zo-z)il-z)-P. (A.7) 

It is important to note that the functions d{m, e), d'{m, e) and constants oq, uq in Table 20 

*The modifications introduced above do not spoil any of the cancellations that take place among the 
original subtraction terms. Hence the modified counterterms are still correct regulators of all kinematic 
singularities. 



[dp£^] « (1 - 
[dpS^^] - (1 - 
[dpg] « (1 - 

and finally 

[dpSr^] = [dpS+i][dpg.]«(l- 



- 22 - 



Doubly unresolved counterterms 


Subtraction term 


Momentum mapping 


Function 


Mo,o) 

irs 


{p} ^ {pjJ;''^' 


/(ofo, a,rs, d(m, e)) 


ir;js 


{p} {p}*;'-^^^' 


/(ofo, Qir + ajs,d{m, e)) 


^ /io(0,0) 




/(ao,ajr,rf(m,e)) 
x/(yo,yjQ,d'(m,e)) 


o(0,0) pa 9(0,0) 

r c("'0) r 

^ ^c: c(o.o) 




f{yo,yrQ +ysQ - yrs,d'{m,e)) 



Table 20. The modified doubly unresolved subtraction terms are obtained from the original coun- 
terterms (first column) by multiplication with an appropriate function (last column). Also shown 
are the momentum mappings used to define the subtraction terms (middle column). 



are the same as those in ah other modified subtraction terms, discussed in Refs. [25, 31]. 

The form of the exponents d{m,e) and d'{m,e) is actuahy fixed by the prescription 
adopted in Ref. [25] (see in particular Eqs. (3.2), (3.12) and (3.13) therein) and the require- 
ment that the modified subtraction terms should still correctly regularise all kinematic 
singularities. In fact, we must have 

d{m, e) = 2m(l - e) - 2do , and d'{m, e) = m(l - e) - d'^ , (A. 8) 

where do and are the same constants which appear in eqs. (3.2), (3.12) and (3.13) of 
Ref. [25] i.e., 

do = Do + die, and do = Do + d[e, (A.9) 
where Do, D'q > 2 are integers, while di, d'l are real. 



B Spin-averaged splitting kernels 

In this Appendix, we list the spin-averaged splitting kernels. Although some of these 
already appeared elsewhere, for the sake of completeness, we give all the splitting functions 
of our computations. 

B.l Two-parton kernels 

The azimuthally averaged two-parton splitting kernels are well known, 

2 



rW(z) = 2C7a 



1 1 

z 1 — z 



2 + z- z^ 



Pil\z;e)=Tn 



1 



1 



(. - z^) 



2 + (l-e)z 



(B.l) 
(B.2) 
(B.3) 
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In our convention the ordering of the labels on the splitting kernels is usually meaningless, 
but in Eq. (B.3) z means the momentum fraction of the second label, i.e., pfq{^z\t) = 
Pq^ iX — z;e). The other two cases are symmetric with respect to z ^ 1 — z. 



B.2 Three-parton kernels 

The spin average of the splitting kernels was computed in Ref. [52], however the forms 
presented there for gluon splittings are not suitable for us, as it was explained in Ref. [22]: 
in the gluon splitting kernels Pg^qrqs and Pgig^gs^ the terms that depend on the transverse 
momenta must always be written in the form /c^ ^ A;^ fc/^J-j '^J-.fc '^^'^ equal to j). Oth- 
erwise the collinear behaviour of the counterterm cannot be matched with that of the singly 
collinear counterterm in the singly unresolved region of phase space. The correct azimuth 
dependence can be achieved by substitutions according to the following replacements, 

(B.4) 

2 j^J-,fc ~^ y^jk ~^ 2 ZjZkSj]f_i — ZiSj(iif_'^ — ZjSi(^jf^^j ■ 



k\ i ■ k\ 



,k 



where {k, 1} = {i, r, s} \ {j} and j can be i, r or s. With these forms azimuthal averaging 
amounts to the simple substitutions, 



(B.5) 



\ k±j ■ k^,k I 2(1 - e) ' 

where (...) denote spin averaging, and j = k \s also allowed. 

For quark splitting into unequal and equal quark flavours we have. 



and 



Pq%qs i{^i,kh Sjk}; e) = 'iPqmWs ii^hkh Sjk}; e) + -Pg'^^g/ {{zj,ku Sjk}; e) , (B.7) 
respectively (hence Eq. (3.9)), where 



— Pq^q'rq's^{^3M^^jk]]^) = CfTrI 



^r.is ~1~ ^s.ir ^rsy^vAs ~1~ ^sAr) ^irs^rs 



+ s (z ■ +z ■ ~ z ■ +z ■ \ ~ '7~ 

+ (r o s)| 

(B 
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and 



2 J 



Sr. 



l + z: 



r,is 



2z, 



1 ^7\iR 



+ l + Zr 



2z, 



+ Zr 



+ Zs 



1 + <is 



e 1 + 2 



Zi rs ~l~ Zr zs 



+ s)} . 



For splitting into a quark and a gluon pair we find, 



where 



(B.9) 



Pi-LsiUjM^ ^Jkh e) = CfP(J),^({z,- H, sjkh e) + C^P^^^^^lUz^^u, s^k]; e) , (B.IO) 



1 ^s^ir ^1 ■^r 

2 

T,is \ Zr is '\~ Zs ir 



+ 



1 - Z 



r,is ■^s,ir 



Sir Si, 



2 + (1 - e)Zs,ir + 1 



+ 



1 - e 



e - —(1 + e)(3 - Zr^is - 2zs^ir) - —(1 - e) 

Sir 5' 



+ (r f-)- s) 



(B.ll) 
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and 



CaCf 



(1-6) 



^ir^s,ir ~^ ^is^r,is ^ir^it. 



1 (1 ^r^is) ~l~ (1 ^s,ir^ 



1 



1 + ^2 



+ 



1 



+ 



1 



(1-6) 



1 



i,rs _j_ ~^ 



ir J 



Zr.7 



^r,is 



^s,ir) 



2 



1 



+ 



+ 



1 



(1 - ef 



1 + (1 — Zs^ir)"^ 4 '^Zs^ir + Z^^^^ Z^^is 


6 


f ' ( 


Zr,is Zf^ig + Zg^ir 


+ 2 





Zr,is ^ Zg^ir 



s,ir 



^rAsi^ Zs,ir) \ Zi 



Zg.ir{.Zr,is ~1~ -^s.zr) 



+ 



Sir Si, 



+ 



+ 



+ {r ^ s)} . 



(B.12) 



We call attention to the normalisation (in colour space) of the abelian and non-abelian 
parts of the splitting function in Eq. (B.IO), which is the same as in Ref. [20]. Notice that a 
factor of Cf, respectively Ca, is made explicit in the definition of the abelian, respectively 
non-abelian, part as compared to the complete splitting function. However, we prefer to 
define [C^y^Jj^j^j^ and [c||?^]j^j^j^ to be dimensionless in colour space, hence the factors of 
Cf and Ca are not made explicit in the definition of these functions, see Eq. (4.2). Then, 



Eq. (4.1) must be interpreted with some care when computing [C-^]]^^^^^^ or [0^^^]^ ^"^^. In 



,(0)n(nab) 



particular, we must remember to include the factors of Cp and Ca explicitly with P 



(ab) 



and Ppf^f^, i-e., we must set Pf^j^f^ 



CpPf p ^ r or Pf \ f 

Jijrjs Jijrjs 



fi fr fs 

CxP'ff^} to obtain the 

Jijrjs 



correct normalisation. 

For gluon splitting into a gluon and a quark pair we have 



PgiLqs iUj,kh Sjk}; e) = CpP^Jg^ {{zj,kh Sjk}; e) + CAPg^qJ, {{zj,ku Sjfc}; e) , (B.13) 



where 



irs 



CfTr 



1 1 
72" 



1 



1 6 Sif -\- Sis 



1 1 



2 Si'ps 

+ Zi 



1 ^ SirgSir^ig 
(1 (^)iZi^rsSrs ^i,rs(l ^i,rs)'S2rs) 



+ (r O s) 



(B.14) 
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and 



CaTr 



,ir 



+ 



+ 



+ 



1 / ^r.is ^r.is 



+ 



1 1 — Zr 



1 1 

+ 



1 1 - e 



2 



1 1 ^Zr.i.^Z. 



r,is'^s,ir 



1 



+ 



1 



1 € Zi^rs\Zr,is ~\~ Zs^ir 
1 



is) Sirs 

2Zr 



+ l-e 



{2Zr,isZs + Zi 

,rs^rs ^r,is^is ^s^ir^i 



2(1-6) 

{2Zr,isZs,irSirs + Zi 



+ 



1 



2(1 -e) 2 



-2^i,rs(l Zi^j.g)Sifs) 



+ r O 



(B.15) 

Note again that the normahsation (in colour space) of the abehan and non-abelian parts 
of the sphtting function in Eq. (B.13) matches that of Ref. [20], i.e., a factor of Cp and Ca 
is made expHcit in the definition of the abelian and non-abelian piece respectively. Hence, 
the comments below Eq. (B.12) apply in this case as well. 
Finally for splitting into three gluons we have 

^irs 

.2 1 f(l-e 



C 



ATI 



+ 



+ 



2'Sij'sZ'p^isZs^ir 
SrsZi^rsi.^ Zi.r 
S 



+ Zs 



+ 



ZrAs H" Zg 



2(1 - Zi^rs) + "^zl^s _ 1 - 2Zi,rs{l " Zj^rs) 



1 - Z 



i,rs 



SirSr 



SirsZi^rsi^ Zi^rs) SrsZi^rs r-. r\ 
Zr^is 7Z ^ v-L ^Zg 



2^s,ir(l -2-5, ir) 



SirsZs,ir(,^ Zs,ir) Si^Zg^ir , 
~r Zy^is ^ yL ^Zi 



(B.16) 



^irs I ''^Zi_^ysZs,ir H" '^Z^^isi^ Zf 



1 1 2^7.^^5(1 Zj- 



(1 ■2-i,rs)(l -2-5, ir) 



^i,rs'^s,ir 



22^5, jr(l Zg^ir) 2^ 2SirsZi^rsZs,ir SifZg^ir S^gZi^rs "t" SigZf^ig 



+ (5 permutations) j . 
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C Integrating the triple collinear counterterm 



C.l Master integrals 

The triple collinear momentum mapping leads to the exact factorisation of phase space as 
given by Eq. (2.12), where the two-particle factorised phase space can be written as 

[<iP2Z\PriPs,Pirs;Q)] = da (1 - a)^^™"^)^^"')"^ d(/>3 (pi , , i P(irs) ) ,^ ^, 

ZlT (C.l) 

X e(a)G(l - a) , 

with p^^,^^^ = (1 — a)pfrg + aQ^. When writing Eq. (4.1), we have used the fact that the 
spin correlations present at the level of the factorisation formula vanish upon azimuthal 
integration via the usual arguments, since k'^ j k Hi k = i, r and s) as defined in Ref. [22] 
is orthogonal to p^^^. Therefore, the integrals of the spin-dependent and spin-averaged 
splitting functions are equal. 

The spin-averaged triple collinear functions pj^j ^ depend on six variables: three two- 
particle invariants: Sj^, Sis, Srs and three momentum fractions: Zi^^s^ Zr^is and Zs^ir, and are 
given in Appendix B. When counting the number of independent kinematical structures 
in Eqs. (B.8)-(B.16), we make two observations. Firstly, of the six variables, only four are 
independent, Zi,rs + Zr^is + Zs^ir = 1- Secondly, from Eq. (C.l) 

it is clear that the factorised phase-space measure [dp2^j^ {pr,Ps,Pirs] Q)] is fully symmetric 
under permutations of the indices {i, r, s}. This permutation symmetry relates integrals of 
different terms, 

/ [dP2Z\Pr,Ps,Pirs;Q)]P{Pi,Pr,Ps) = / [dp2Z\Pr ^ Ps, Pirs] Q)]P{Pa{i) , Pa{r) , Pa{s)) , 
J2 J2 

(C.2) 

where cr G 53, with P an arbitrary function of momenta (a term in the splitting functions). 
This further reduces the number of independent structures. 

Making use of the constraints among variables and the permutation symmetry, we 
find 46 structures to integrate^, that can be grouped into five classes, 



p{j,k,l,m 

p{j,k,l,m 
T-,(j,k,Lm 



Pr 



fj fk A m 

A J- . z"^- 



i4<„,(l-z.,,.)-, (C.3) 

(\ - z ■ Vz"" 



5 ~ ^is^rsi^ Zr^is) (1 + ^s.ir) : 



where we introduced the scaled two-particle invariants tki = Ski/sirs {k, I = i, r, or s). 
The exponents take values as given in Tables 1-5. With these five classes, we can give 
new forms of the spin-averaged splitting functions that lead to the same integrated triple 



This basic set of integrals is not unique, nor do we claim that they are linearly independent. 
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collinear subtraction terms, 



^ f'^/ n=l j,k,l,m 



where the constants CLj.f^f^ are given in Eq. (4.4), while Cj'^'''j^'-'™ in Tables 1-5. Thus we 
are left with five types of integrals, 

X f{ao,a,d{m,e)) , (n = l,...,5). 

C.2 Explicit representations 

In order to write the integrals in Eq. (C.5) explicitly, we must choose a specific parametri- 
sation of the factorised phase-space measure. To begin, we choose the scaled two-particle 
invariants t^i = Ski/sirs, {k,l = i,r or s), and Vr, Vg, defined by 

_ Zr^is - zj ^ _ Zs^ir - zj ^ 

- s+) _ S-) ' - J+) _ J-) ' ^^-^^ 

as integration variables. The momentum fractions take values between 

/ ,\ a + (1 — a^x^— I ^ a 

zi-'^ = {l-Us)-^ 4-) = (i-t,,)^ ^ (C.7) 

la + (\ — a)x~ 2a + (1 — ajx— 

^ ^ irs ^ ' irs 

and similarly for z^^^ with r -H- s. Clearly, Vr and Vg are simply momentum fractions 
rescaled to take values between zero and one. In terms of these variables, we have 

a + (l — a)x~Vk 
Zk,ii = 1 - til) ^ ^ ' , kj = r,s. C.8 

2a + {l-a)x^ 

Using the variables t^i and Vk, the factorised phase-space measure reads 

X yl'J^'^dtir dtis dtrs dvr dvs dyirs6{a{a + (l - a)x~) - yirs) 

X 5(1 - tir - Us - trs){l - Ur){l " iis)[(4t^ " trs){trs " 47^)]"^"' 

X 6(1 - tir)&{tir)e{i - Us)e{tis)e{tit'^ - t..)e(t.. - 4;)) 

X 9(1 - t;,)e(z;,)e(l - z;,)e(i;,) , 

(C.9) 

where 

= (l-t,,)(l-t,,)TW with = [Vvril - Vs) ± - ^r)] ^ • (C.IO) 

Eq. (C.9) is not yet in a very useful form for computing the integrals because it contains 
constraints in the form of nontrivial functions. In order to isolate the singular behaviour 



-29- 



we map the region of integration onto the unit hypercube such that physical singularities 
are on the borders. Solving the constraints in a particular way, we find a parametrisation 
over the unit hypercube where all physical singularities are on the border except Zs = 0, 
which introduces a line singularity (see below). We find 

X dtis dTrs dvr dws dyirs 6{a{a + (1 - a)x~) - yirs) 

X [Ws{l - Ws)]~^~^ (1 - Trs + TrsUs)"^'^^' 
X 6(1 - tis) eitis) 9(1 - Trs) Q{rrs) 
X G(l - Vr) e{Vr) 9(1 - Ws) e{Ws) ■ 



In terms of the variables tis , Trs , Vr and Ws we have 



and 



a + {1 — a)x^Vr 

^r,is — (1 tis}' 



2a + (1 — a)x'^ 



i Tys \ TrsT^is 



(C.ll) 



, _ (1 '7Vs)(l tjs) , _ Trstis{l Us) , „ . 

-L Trs ~r Trs*^is -L Trs ~r TrsT^is 



a + (1 - Q;)x-;j[Trs(l - Vr) + Ur(l " T^s) - 2^Trs{l - Trs)Vr{l - Vr)^^ " 2Ws)] 

2a + (1 — a)x^ — 

(C.13) 



We see that = corresponds to a = 0, = 0, 1 and Trs = Vr, hence the line singularity. 
Since the only integral involving 1/zs is X2c,4, this is the only place where the line singularity 
has to be resolved. 

Using the parametrisation of Eq. (C.ll) and Eqs. (C.12) and (C.13), we find the fol- 
lowing explicit parametric integral representations of the basic integrals. 



X^^^y''^\x~,e;ao,do) = 2 ^^{1- 2e) ^ ^ dUsdTrsdvrdws 



JO 



X a 



-l-2e 



(1 _ a)2'^o-3+2. + a)x~)-i-2^ (2a + (1 - a)x~)- 



X i (1 - t..)^+'+™+^-'^ Ttr (1 - T..)^-^ (1 - Vr)-' 

X u;;^"' (1 - Ws)-^^-' (1 - Trs + r,, t..)-^'-'-2+2^ (a + (1 - a)x~Vry 



(C.14) 
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j-ij^k^ ^'"'^(a;~,e;ao,'^o) = 2 ^y(1- 2e) / "^"Z ^^^^ ^'^^'^ "^^^ "^"^^ 
X (1 - a)2^o-3+2^ (a + (1 - a)x~)-'~'' (2« + (1 " 

X (1 - T-/ (1 - TrsV-^ (1 " VrT^ 



"is \^ "ts; 'rs \^ "'sj 

X W^"' (1 - Ws)-'^-' (1 - Trs + Trs («+(!- 

X |a + (1 - a)x~[(l - Trs)Vr + Trs{l - Vr) - 2^/ Trs{l - Trs)Vr{l - Vr){l - 2Ws)]^ , 

(C.15) 

2-0,M,m)^^_ , e; ao, c^o) = 2~^' Jr(l^- 2e) / J ^^""'^ ^^"^ 
X (1 - a)2'^»-3+2^ (a + (1 - q)x~)-^-2^ (2a + (1 - 0)2;—)-™-^ 

X i ■ (1 - r-f (1 - T,,)^'-^ (1 - Vr)-' 

X W^"' (1 - Ws)-^-' (1 - Trs + TrsUs)-'-^-^+^' 

X [(1 + tis)a + (1 - a)x~(l - (1 - tis)vr)r 

X |a + (1 - a)x^[(l - T,.s)t;r- + Trsil - Vr) - 2y^Trs{l - Trs)Vr{l - t'r)(l " 2Ws)]^ , 

(C.16) 

^2C,4~^ '"^^(^i^^ ' «o, rfo) = 2"^^ _ 2e) / / "^^^^ '^^'■'^ "^^'^ ^^"^ 
X (1 _ afdo-3+2e + a)x~)-i-2^ (2a + (1 - a)x~)2 

X 4-2^ (1 - r,^r (1 - rrsY-' vr (1 - Vr)-' 

X W^"' (1 - Ws)-'^~' (1 - r,, + Trstis)-^-''-'^^' 

X [(1 + ti,)a + (1 - a)x~(l - (1 - tis)?;,-)]"^ 

X |a + (1 - a)x~[{l - Trs)Vr + Trsi^ - Vr) - 2y/Trs{l - Trs)Vr{l - Vr){l - 2Ws)]j 

(C.17) 

and 

^2C^5~^~^~^^(^l?^^^'i"o,do) = 2^"^" —x:r^ I da I dtisdrrsdvrdw, 



7rr(l - 2e) 



JO 



X _ ^)2do-3+2. + a)x^)"i-2^ (2a + (1 - a)x~) 

1 



X (1 - r^, + r^,ti,)"^+2^ [(1 + tis)a + (1 - a)x~(l - (1 - tis)vr)]-^ 



a + (1 - a)x^[(l - Trs)Vr + Trs{l - Vr) - 2-^/Trs{l - rr.s)Ur(l " Vr){l - 2Ws)] ^"1 

^ 2a + (l-a)x~ / 

(C.18) 

The expressions for X2c,n., ^^ = 1, 2 and 3 in Eqs. (C.14)-(C.16) are suitable for evalu- 
ation with general purpose sector decomposition codes as they stand. In X2c,5, Eq. (C.18), 
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the terms in the braces require some care. Although this factor is just 1/(1 + Zs), which is 
finite and hence in principle it can be simply carried through sector decomposition as it is, 
nevertheless depending on the precise internal implementation, it may lead to undefined 
arithmetic expressions in the (symbolic) computation. We find that no such trouble arises 
either with our private implementation or with SecDec, if we gather the two terms in the 
braces over a common denominator. 

2^2C,4 requires special attention, as could be anticipated by the presence of the factor 
l/zs- There are two problems: firstly, the presence of the square roots in the vanishing 
denominator in the braces on the last line of Eq. (C.17) interferes with the treatment of 
overlapping singularities. (This could be solved by simply making the expression square- 
free.) Secondly, as we have indicated, there is a line singularity inside the integration region 
at Trs = Vr- We address both issues by deriving an alternative representation for X2C,4; 
which is free of both square roots and line singularities. The price to pay is that the new 
representation is quite cumbersome. 



r2(l-e) 
"^'r(l-2e)^' 



da 



dt,o dTrQ ds du a 



4-2e .fc-2e 



-l-2e 



X (1 - a)2'^°-3+2^(a + (1 - a)x~)-^"2.(2a + (i _ a)x~)'-H 



irs 



l-2e 



2(1 + tis)a{a + (1 - a)x~)(l + u) 

+ (1 - a)^(l - Trs + TrsU)x~ + tis(l - a)'^(Trs + U - TrsU)x^. 

1 + tis)^a{a + (1 - a)x~) + ti,(l - afx^-^'j (2a + (1 - a)x~;)\l - sfu^ 
+ \a{a + (1 - a)xr^)s{l - uf + (2a + (1 - Ci)x~fu 
X s [2(1 + tis)a{a + (1 - a)x~)(l + u) 

+ (1 - a)^(l - Trs + Trsn)x~ + tjs(l - a)^(Trs + n - rr.sn)x~ I 

(C.19) 

The derivation of this alternate form relies on rewriting the integration over Vr and Wg, 



'■Vr,Ws 



dVrdWsV^'^{l — Vr) '^Ws^ '^{^ — Ws) 2 ^ 

X [a(l + tis) + (1 - a)(l - (1 - tis)vr)x~:;] 

X |a + (1 - a)x-^[(l - Trs)Vr + Trs{l - Vr) 

- 2^Trs{l - Trs)Vr{l - Vr){l - 2Ws)]^ , 



(C.20) 



as an angular integral, 

= 2^+4^(1 + Us)-H2a + (1 - a)x~)-2^ J dO,_i(g)^^-^ 



q){P2 ■ q) 



(C.21) 
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where in a suitable frame, 
1,... 



Pi 



{l-tis){l - a)x~:; 
(l + t,,)(2a + (l-a)x~) 



P2 



[1 — a)x~ 



(the 



2a + (1 - a)x-r^ 
denotes vanishing components) and 



■smx, 



(1 — a)x' 



2a + {1 — a)x- 



■cosx 



= (1, . .'angles'.., sin '!?sin(^, sin t?cos(^, COST?) 



(C.22) 
(C.23) 

(C.24) 



(Here ..'angles'., represents those angular variables that are trivial to integrate, since the 
integrand does not depend on them.) Clearly both pi and p2 are massive and time-like, pi, 
P2 > 0, and pi • P2 > 0. The angular integral can be written in terms of a Mellin-Barnes 
representation according to Eq. (60) of Ref. [56] . Finally, we can perform the Mellin-Barnes 
integrations at the expense of reintroducing two real integrations over the unit interval. 
We obtain a new two-dimensional real integral representation. 



l-2e 



2i+^%e(2a + (l-a)x~)2-2^ 
1 

ds du s"2'n~2' 2(1 + tis)a{a + (1 - a)x^)(l + u) 

+ {1- a)^(l - Trs + TrsU)x~ + tis{l - af{Trs + U- TrsU)x'~ 

(1 + tisfa{a + (1 - a)x~) + ^^,(1 - afx^^'^ (2a + (1 - a)x~)^(l - s) V 
+ [a{a + (1 - a)x~)s{l - uf + (2a + (1 - a)x~r^fu 
X s [2(1 + tis)a{a + (1 - a)x~r )(1 + u) 



-I- (1 - afil - Trs + TrsU)x~ + tis(l - a)'^ (Trs + U- TrsU)x~ 
\ / \ / ^rs \ / \ / ^rs 

free of square roots and singularities inside the integration region. 



-l+e 



(C.25) 



D Integrating the double collinear counterterm 
D.l Master integrals 

The double collinear momentum mapping leads to the exact factorisation of phase space 
as given by Eq. (2.13), where the two-particle factorised phase space can be written as, 

[dptrii'^(j>r,Ps,Pir,Pjs; Q)] = da d/3 (1 - a - /3)2(— i)(i-^)G(a) e(/3) 6(1 - a - /3) 

s- s-TjQ (D.l) 

X -^^d-Cl)2{pi,Pr;P{ir)) ^^#2 (Pj , ; ) , 

where p^.^^ = (1 - a - /3)pf^, + aQi^ and p^^.^^ = (1 - a - /3)p^^ + 

When writing Eq. (4.5), we have used that spin correlations generally present at the 
level of factorisation formulae cancel after azimuthal integration, since ■ and /c'f • as 
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defined in Ref. [22] are orthogonal to and p^^ respecitvely. Examining the actual form 
of the spin-averaged Altarelli-Parisi functions, and using the symmetry of the phase-space 
measure under the momentum exchanges pi ^ Pr, or pj -H- ps, and under exchange of the 
pairs {pi,Pr) ^ (Pj^Ps), we see that we must compute the following integrals, 

^2c'6(^i;^'^7;i^'«0,'^0) = (^"-P-Q^'^ [ [dP2!r^''\Pr,Ps,Pir,Pjs;Q)] — — 



(D.2) 

^r,i4j/(«o, a + /?, d{m, e)) , 
where k, I = —1, 0, 1 and 2. By symmetry, we need only consider the cases when e.g., k < I. 

D.2 Explicit representation 

To write the integral in Eq. (D.2) explicitly, we must choose a specific representation of the 
factorised phase-space measure. We choose the scaled two-particle invariants yir and yjs 
and rescaled momentum fractions Vr^i and Vgj, 



^sj = -TT. TZS , (D.3) 



as integration variables. Here 



a + (1 — a — /3)x- 



a 



2a + (1 - a - /3)x~ ' ' 2a + (1 - a - f3)x- ' ^ ^ 

with similar expressions for Zs^^ and zi ^ with a — )■ /3 and ir — )■ js. In terms of a, /3, Vr,i 
and Vsj- we have 

yir = a{a + (1 - a - , yjs = l3{a + (1 - a - /3)xj'^ ) , (D.5) 

_ a + (1 - a - /3)x~:vr,i _ /3 + (1 - q - /3)x-7;,j 



2a + (1 - a - ' 2/3 + (1 - a - /3)x~ 

and for the factorised phase-space measure we obtain 

[dP2■m'\Pr,Ps,P^r,Pjs;Q)]= ((^<3~^') (Q^)^ Q ^TsQ 

X dad/3 (1 - a - /3)2(™-i)(i~^)e(a)G(^)e(l - a - /3) 

X dyird?;.r,j[yir'fr,i(l - Vr,i)]~^6(^a{a + (1 - a - /3)x~) - y^,.) 

X dyjsdvs,j[yjsVs,j{'i- - Vs,j)]~^s(^l3{/3 + (1 - a - /3)x~) - yjs 

X G(Wr.i)0(l - Vr,i)Q{Vs,j)Q{l - Vsj) ■ 



(D.6) 



(D.7) 
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Using the parametrisation of Eq. (D.7) and the expressions in Eqs. (D.5) and (D.6), we 
find the following parametric integral representation for the master integral, 

^2C,6(^i;^'^J5'^'"o,4) = a^i^ra;— j dadfS Q{ao - a - (3) j dvdu 

X a-l-^ /3-l-^ (1 - a - /3)2'^0-2{l-.) _ ^ye ^-e (-^ _ ^^-6 

X [a + (1 - a - /3)x~]-'-' [/3 + (1 - a - ^)x~]-'-' (D'^) 



X 



a + {1 - a - I3)x~v\'' / P + {1 - a - P)x~uV 



2a + {l-a- P)x~J y2P + {l-a- P)x~ 



Eq. (D.8) is directly suitable for treatment with general purpose sector decomposition codes. 

E Integrating the soft collinear counterterms 
E.l Master integrals 

The consecutive collinear and soft mappings lead to the exact factorisation of phase space 
as given by Eq. (2.14), where the one-particle factorised phase spaces can be written in the 
following form. For the collinear mapping we have 

[dPl-l+M^P^r-^ Q)^ = - «)'"^(^-^)-^ ^ dMP^,Pr■,Pi^r)) e(a)G(l - «) , (E.l) 

where p'^-^^ = (1 — a)p^j. + aQ'^. For the soft mapping we find 

[dpStfe, K; Q)] = dy(l - y)(™-i)(i-^) ^ d^Ps, K; Q) 9(^)6(1 - y) , (E.2) 

where y = ygg and the momentum K is massive with = (1 — y)Q'^. As the notation 
above indicates, a and y are integration variables. 

To integrate Eqs. (4.8)-(4.13) over the factorised phase spaces of Eqs. (E.l) and (E.2), 
we can pass to the azimuthally averaged Altarelli-Parisi splitting functions, since the az- 
imuthal correlations generally present in Eqs. (4.8)^(4.13) vanish after integration by the 
usual argument (the transverse momentum k^^^ in the splitting kernels is defined to be 
orthogonal to the parent momentum p^^). The Altarelli-Parisi functions can be expressed 
as linear combinations of powers of momentum fractions, so we need to compute integrals 
of the form. 



(47r) 2e\ /"r- {ir) , - .r^\^\A /OM C I \ 2 1 - Zg^ir 2 Z. 



o Q ] mZ+M,p^r■, Q)] mLiPs; Q)] ^s.kis) , 

X — 4,j /(«o, a, d(m, e))f{yo,y, d{m, e)) . 



We perform the integration over [dp[*^_^]^] first, so we need to write the integrands in 
terms of ps and tilded momenta. The singly collinear mapping of the momenta rescales 
the invariants which do not involve the collinear pair, 

Vjk = (1 - afVjk yjQ = (1 - a)y.Q , (E.4) 



(E.3) 
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while those which involve the collinear momenta become 

y{ir)k = (!-")[(!- «)y i^Tfc + " y-kq] y^g = (i - ")y I^q + 2a • (E.5) 

The successive singly soft mapping of the momenta rescales those two-particle invariants 
that do not involve the soft par ton, 

y-jk = " y^yfk ' (^•^) 

while those involving the soft parton s become y--- = Finally, we have 

y-jQ = (1 - y)yjQ + yjs > (e-7) 

for j 7^ s, while ysQ = y is simply an integration variable. Therefore, the integrands are 
expressed as follows. For the eikonal factor, we find 



^Ms) = j^:^s. (E.8) 



if j and k are distinct from (ir), while 



2 ^ (1 - a)2 [(1 - + as,Q]s-,^ ^ (1 - a)^ [(1 - a)s~. + as.g].^, ' ^ " ^ 



if e.g., j coincides with (ir). Finally, we have 

2 1-Zs,ir_ 2 s^ir)Q _ 2 2a + (1 - a) [(1 - y)^-:^ + . 



S{ir)s Zs,ir s^Q (1 " a)^ [{1 - a)s~ . + assqjss 

and 

Sjs Zs,j Sjs SsQ (1 - SsQ 

Hence, we define five soft collinear master integrals. 



(E.IO) 



(E.ll) 



Q'^) JWl■l+liPr,P^r■, Q)] [d^^i^fe; Q)] 

1 1 - y 1 I 

2 (1 — aj^ J Sjr ' 

X /(ao, a, e))/(2/o, c^lm-, e)) , 



(E.12) 



^ 1-y (l-a)s~fe + agfcQ 1 
(1 - q)2 [(1 - a)s~. + aSsQ]sy Sir 
X /(ao, a, (i(m, e))f{yo, y, d{m, e)) , 



(E.13) 
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^2(S,3(^i;^; ^' "0, yo, do, do) 
'(47r 



^g2e^ ^[dpg+i(p.,K.; g)] [dri^i^fe; g)] 



(1 - q)2 [(1 - a)s~. + as5Q]s^,. Sir 
X /(ao, a, d{m, e))f{yo,y, d{m, e)) , 

^2 \ 2 



(E.14) 



{^Q'') JWi^i+M-^P^r-, Q)] [dritfe; g)] 

2 2a + (l-a)[(l-y)s~Q + .~,] 1 ^ 



(E.15) 



(1 - a)2 [(1 - a)s~^- + aSs-Q]s^ 
X /(qo, a, (i(m, e))f{yo,y, d{m, e)) 



and 



4c?,5(^7;^' ^7;^ J,q; «o, yo, do, d'o. 



^Q'^) /Jdpi;:^+i(p.,p„.; Q)] mUPs; Q)] 



(E.16) 



(1 - aYs--~ SsQ Sir 

X /(ao, a, d(m, e))f{yo, y, d{m, e)) . 
We need to compute these integrals for / = —1, 0, 1, 2. 
E.2 Explicit representations 

Eqs. (E.12)-(E.16) are very similar to the iterated collinear - soft collinear integrals we 
computed in Ref. [31], and we employ the techniques of Ref. [31] for computing the integrals 
-^2(5 n- begin by recalling the specific representation of the factorised phase-space 
measures in Eqs. (E.l) and (E.2). For the singly collinear measure in Eq. (E.l), we choose 
the scaled two-particle invariant yir and the rescaled momentum fraction Vri, 



(E.17) 



as integration variables. Here 



2a + (1-0)2;— 2a + (l-a)x-^ ^ ' 



In terms of a and Vyi-, we have 



a + (1 — a)x-^Vri 

a(a + (1 — a)x^) , and Zri = ; ^ — (E.19) 

^ ^ ^ ' 2a + (l-a)x- ' ^ ^ 
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while the phase-space measure [dp^^!2,+iiPr,Pir',Q)] reads 

X 5(^a{a + (1 - a)x-r:^) - yij^Q{a)Q{l - a)Q{vr,i)Q{l - Vr,i) ■ 

(E.20) 

Turning to [dp^'J^fe, K; Q)] in Eq. (E.2), we use the (scaled) energy and the angular vari- 
ables of ps in some specific Lorentz frame to write the two-particle phase space d4>2{ps,K; Q) 
explicitly. In particular, we work in the rest frame of Q^, where = • • •) and ps is 

parametrised as, 

= Es{l, ..'angles'..,sinT?sin99sin77, sin?? sin (/3 COST/, sin "i? cos cos "i?) . (E.21) 

The specific orientation of this frame is chosen below according to what is most convenient 
for computing each integral. However, independent of orientation, in terms of the scaled 
energy-like variable, 

_ 2ps.Q _ 2Es 

- -w~ ~ W ' ^ ^ ^ 

and the angular variables 99 and 7], the two-particle phase space d(p2{Ps, K;Q) reads 

dMPs, K; Q) = ^g^5,(-22^6)de, eY^'5{y - e,) 

X d(cos -d) d(cos 99) d(cos 77) (sin ■d)'"^" (sin (/j)"^"^" (sin 77)"^"^^ . 
We are now ready to display the master integrals. 

Integrated soft coUinear counterterm for i, j, k distinct. This case leads to the 
integral X2cs,i- We choose the orientation of the frame such that the momenta appearing 
in the integrand take the following forms, 

# = £^j(l, . . . , 1) , p^ = £^fc(l, . . . ,sinxfc,cosxfc) ) ,^ 

. (E.24) 
Pir = Eir{l, . . . , sin cl)ir sin Xir , cos (pir sin Xir , cos Xir ) • 

Then, expressing the various invariants in terms of the integration variables in the chosen 
frame of Eq. (E.24), we find 

yjcs = ^V'kQyi^ ~ Xk sin ■!? cos V9 - cos Xk COS 1?) , (E.26) 

yirs ~ 2 ^ ir ~ "i^*^ ^™ ^ ^™ ^ ^ ~ ^™ ^ ~ ^) ' 

(E.27) 

where 

cosxfc = cosx(>jfc^Q) , cos Xir = cosx(^j~ q) 

Y-rr ^+Y-r~ ^-Yrr- ^-2Y-rr ^Y-r^ ^ (E.28) 



COS ( 
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with 

cos xiY) = I - 2Y , smx{Y) = 2^Y{l-Y). (E.29) 

The factor y^. in the denominator of the eikonal factor vanishes at cos ip = I and 
COST? = cosxfc- Hence, the integrand has a hne singularity that we remove by performing 
partial fractioning as in Refs. [31, 57], 



s-js^ks ^'"'^y^y-js ^yksJ\yjQ y-kQ 

Then we find 



yo 

dyyi-2^(l-y)'^o-i+^ 



"1 

X / d(cosT?) d(cos(/9) d(cosr/) (sint?)"^^ (siny?)"^"^^ (sin?/)~^~^^ 







-1 

X / daa-i-^(l-a)2'^"-3{a + (l-a)[(l-y)x~+y~,]}-^-^ 



(E.31) 



X / dtif ^(1 — v)' 







g + (1 - a)[(l - y)xj^ + y-^^v 
2a + {1 - a)[{l - y)x~ + y 



X [(1 - y)x~ + J hr^ + TT-^ + 



1 



^ 2yjg- 2yy y \^ y--Q y^Q J 

where y^-, yj,~ and are understood to be functions of the integration variables as given 
in Eqs. (E.25)~(E.27). We also used y^q = x^^, and performed the integration over with 
the (5-function in Eq. (E.23). The choice of frame in Eq. (E.24) is convenient for integrating 
the first term in the partial fraction, while the second term is more straightforward to 
integrate in a frame where k and j are exchanged. Upon performing this exchange, we 
find that the functional form of the master integral is unchanged, and we must simply 
interchange Y~~ „ with Yr~ „ to obtain one term from the other. 

The integral 1^2CS i given in its most general form in Eq. (E.31) first appears in 
electron-positron annihilation into four or more jets at NNLO. For a three-jet computation, 
the r/ integral is trivial, since momentum conservation forces the three final state momenta 
to be coplanar, and hence we find (j)ir = or vr in Eq. (E.24). We display the resulting 
simplifications. Choosing = vr, Eq. (E.24) becomes 



= £^j(l, . . . , 1) , = £^fc(l, . . . ,sinxfc,cosxfc, , 

Pir = Eir{l, . . . , - Sm Xir, COS Xi 



rr J 1 
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which leads to 

y'js = \yjQyO- - cos ^) , (E.33) 

y'ks = ^2/fcQy(l - smxfcsini9cosv3 - cosxfccosi?) , (E.34) 

yi^s = \^l^y(^ + Xir sin cos - cos Xir cos t?) , (E.35) 

with cosxfc = cosx(^~^ g) and cosxir = cosx(^~~ g) as in Eq. (E.28). Then the inte- 
gration over T] can be performed using 

^ d(cos7?) (sinr?)-2-2^ = . (E.36) 

Due to momentum conservation, only two out of the five kinematic variables in Eqs. (E.33)- 
(E.35) are independent. E.g., in terms of y~r „ and „, we have 



jir,Q'^\ j ir ,Q ' \l j k jk 



3k,Q\J ]ir,Q^ 3tr,Q> 3 ir ,Q\I 3 k 3k,Q' 

(E.37) 



3 ir 3 ir ,Q' 

y'ko = H = 1 ; , (E.38) 

3k,Q\l 3 %r 3 '^'r ,Q' 3tr,Qy 3 k ,Q'^ 3k,Q' 



^3k,Q^^ ^jk,C 



3k,Q\J 3ir,Q^ 3^r,Q^ 3 ,Qy 3 k ,Q'- 3 k 

Note that Yj-^ ^ is also easily expressed with Y-j ^ and Y~~ 

k ir ,Q 3 k ,Q 3 ir 



(E.39) 



k ir ,Q 3 k ,Q 3 ,Q 3 k ,Q 3 tr ,Q V 3 k 3k,Q' 3 tr 3 ,Q' 

(E.40) 

The physical region in y~r „ and Y~~ „ is simply given by < Y~r ^,Y~~ „ < 1 and 
^~3k,Q+^~jl^,Q > ^• 

In the case of two-jet production, this integral does not appear at all because we do 
not have three independent hard final state momenta. 

Integrated soft coUinear counterterm for j = {ir). This case leads to the integrals 
1-2CS 2 ^^'^ '^2CS 3 • -^^r 1-2CS 2 integrand again has a line singularity, which we can remove 
via partial fractioning. Using 



[(1 - + aSsQ\sj^~^ Q2 y2 

V / \ -1 (E.41) 



+ -^\\ 2a + {l-a)x~^^ + 



2a + (l-a)x~^ ^y-ks I \ y^Tr y^kQ 
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we find the following parametric integral representation, 



■p2 /'I ^\ fPo 



X / d(cos'i?) d(cos(/7) (sin-!?)"^^ (sin^j)"^"^^ [(1 - y)a;~ + y 



X 



/ daa-'-'{l - a)2'^"-3{a + (1 - «)[(! - y)x~ + y-.]}' 



l-e 



1 , yy-kQ 



X / dt;t;"^(l - v)~^\a + (1 - a)x~Y~r J ^ h 

2a + (1 - ajx^ — ^ "^Vks 

[ '^'y^rf yy-kQ) V2« + (i-a)[(i-y)x~ + y~,], 

(E.42) 

The two terms in the partial fraction are evaluated most conveniently in two different 
frames. In the first one 

Pfr = -^ir-(l,---,l), = ^fc(l,...,sinXfc,COSXfc), (E.43) 

and hence 

= -x—wfl - COST?) , (E.44) 

y-ks = \ykQyi'^ ~ X/c sin •ecosip- cos cos ^9) , (E.45) 

with cosxfe = cosx(^~^ q)- In this frame, the y^. = singularity is at the border of 
integration. Instead, in the second frame, the y^~ = singularity is at the border, 

= . . . ,1) , p'^^ = Eir{l,...,SmXir,COSXir), (E.46) 

(note that clearly cosxir = cosxfc) and hence 

yfcs = ^2/fcQy(l-cos^), (E.47) 
yi:fs = \^l^y(^ ~ sinxirSint?cos(^ - cos Xir cost?) , (E.48) 

with COSXir = cosx(>~fc q). 

As for 1265 3' most conveniently evaluated in the frame defined by Eq. (E.43), 
where we can integrate over ip using 

r2(l-e 



27rr(l - 2e) J_i 



J d(cos if) (sin ip)'^-^' = 2-^'^' . (E.49) 
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Then we obtain the following explicit integral representation, 



^£,3(^7^; «o, yo, do, = -2^^ j^" dy y"^'{l - yfo-^+^ j'^ d(cos ^) (sini?)" 



2e 



/•ao 

X [(1 - y)x~ + y~^] / daa-^(l - af^'-^a + (1 - a)[(l - y)x~ + y~s]}~'~' 
Jo 

f\ 1 fa + {l-a)[{l-y)x~: + y~,]v\' 

X dvv [1 — V) ' 







2a + (1 - a)x~^ V2a + (1 - «)[(! - y):^;~r + y~ 

(E.50) 



Integrated triple collinear — soft coUinear counterterm. For the integral I2GS 4 
use the frame where 

pi; = ^i,(i,...,i), (E.51) 

and hence 

y-^. = -x-r-y(l - cos i?) . (E.52) 
Then we obtain the foUowing exphcit representation, 

Jo 

X y"^d(cos^?) (sin^?)-2^[(l-y)x~ 



X 



/ da a-i-^(l - a)2'^"-3{a + (1 - a)[(l - + y-.]}' 

Jo 



1-e 



(E.53) 



X / dvv [l-v) ^ — 

2a + (l-a)x~^^ 

^ ' ir yx~ 

/a-|-(,i — «;[(,! — y jx— -t- y~ ~\u \ ' 

X 



g + (1 - a)[(l - + 
2a + (l-a)[(l-y)x~+y~.] 

Integrated double collinear soft collinear counterterm. Finally the integral 
-^2CS 5 iiiost conveniently evaluated in the frame where 

= ,1), pf^ = . . . ,sinxir,cosxir) , (E.54) 



and hence 



V~ys = \yjQy(^ - COS ^) , (E.55) 
yifs = \^'i^y(^ ~ sill COS 93 - cos Xir cos t?) , (E.56) 
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with cos Xir = cosx(^^j q)- The exphcit expression for the master integral reads 

■p2 /i -\ pyo 



X y ^d(cos??)d(cos(^)(sini9) ^^(siiK^) ^ "^^ [{I - y)x~ + y~ 

f-ao 

X / daa-i-^(l-a)2'^«-3{a + (l-a)[(l-y)x~+y~J}-i-^ 
Jo 

f\ -en .-eyy-jQl, y]s\ {a + {l-a)[{l-y)x~+y~.^v\^ 

X dvv [1- v) \ l-y + 







(E.57) 



The parametric forms for the integrals 1-2cs n presented in this section are directly 
suitable for evaluation using sector decomposition, after introducing the new variables 

1 — cos ■!? 1 — cos if 

zi = ^ and Z2 = ^ . (E.58) 
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